Nearest-neighbor detection of atoms in a ID optical lattice by fluorescence imaging 



o 
o 

(N 

CD 



Oh' 
I ■ 

-(— > • 



> 

oo 
en 

o 

00 
O 



% 



M. Karski,| L. Forster, J. M. Choi, W. Alt, A. Widera, and D. Meschede 
Institut fiir Angewandte Physik der Universitdt Bonn, Wegelerstrasse 8, 53115 Bonn, Germany 

(Dated: February 6, 2009) 

We overcome the diffraction hmit in fluorescence imaging of neutral atoms in a sparsely filled 
one-dimensional optical lattice. At a periodicity of 433 nm, we reliably infer the separation of two 
atoms down to nearest neighbors. We observe light induced losses of atoms occupying the same 
lattice site, while for atoms in adjacent lattice sites, no losses due to light induced interactions 
occur. Our method points towards characterization of correlated quantum states in optical lattice 
systems with filling factors of up to one atom per lattice site. 
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Neutral atoms in optical lattices have been shown to 
be an ideal system for engineering novel types of strongly 
correlated quantum states. Quantum correlations be- 
tween different lattice sites could be induced with Bose- 
Einstein condensates by precisely adjusting the relevant 
energy scales through controlling the lattice potential 
[y, y . Detection of these novel states was initially only in- 
direct by observing the collapse and revival of the global 
matter wave interference pattern in time of flight mea- 
surements. In contrast, quantum state tomography, as 
well as many theoretical proposals to employ these corre- 
lations for quantum information processing, require sin- 
gle site detection |J, a technically challenging goal for 
site separations in the optical wavelength domain. In 
a different regime, where potential wells are separated 
by several micrometers, single atoms could be resolved 
P, 0. However, in this regime the relevant energy scales 
are not well amenable to control via the external poten- 
tial, therefore the "standard route" for the preparation of 
correlated quantum states sketched above seems to favor 
site separations in the optical wavelength regime. Re- 
cently, single site detection has been reported in such a 
system using focused electron beams from an ultra-high 
vacuum compatible electron gun |_[]. This technique is 
not easily integrated with many current quantum gas 
experiments, in which, in contrast, optical imaging by 
fluorescence light is widely established. The latter has 
seen great success in other fields, e.g., imaging of single 
molecules |7|] . Comparable success with neutral atoms in 
optical lattices, however, could not be achieved to date. 

In this work, we demonstrate the detection of atom 
pair separations down to nearest neighbors in a onc- 
dimcnsional (ID) lattice with optical wavelength period- 
icity. We overcome the previous restrictions imposed by 
the diffraction limit M with a markedly improved data 
quality and reduced noise, together with advanced nu- 
merical processing of fluorescence images. Such a new de- 
gree of precision in detection allows us to directly observe 
light induced atom losses and to distinguish between on- 
site and nearest-neighbor contributions. In contrast, de- 
tecting such loss processes have so far relied on ensemble 
averages in optical lattice systems M, while interacting 



atoms in a type of atom blockade effect have been investi- 
gated in systems where only a single running wave optical 
trap was present |1(|] . Our approach, however, combines 
single atom resolution with single site detection. 

In our experiment, wc load a small number of Caesium 
atoms in the periodic dipolc potential of a standing-wave 
laser field - a ID optical lattice - from a magneto-optical 
trap (MOT). The lattice is formed by two countcrprop- 
agating linearly polarized laser beams (A — 865.9 nm) 
with a waist of wq = 20 /xm and a typical power of 
P = 100 mW providing a trapping potential with a depth 
of {//fee = 0.4 mK, for which atom tunneling is negligible. 




Figure 1: Detail of the experimental setup. Two counter- 
propagating laser beams (a) form the ID lattice. Beam tubes 
(b) shield the objective (c) from stray light of molasses beams 
(d) off the glass cell (e). Apertures (f) strongly suppress the 
remaining stray light. A narrow-band optical filter in front of 
the EMCCD (g) filters the stray light from the optical lattice. 

We illuminate the atoms with a red detuned three- 
dimensional optical molasses at 852 nm which also pro- 
vides continuous Doppler cooling and counteracts heat- 
ing of the atoms. The fiuorescence light is collected by a 
diffraction limited microscope objective (NA = 0.29) Q 
and is imaged onto an electron multiplying CCD (EM- 



CCD) camera, where stray light from the molasses beams 
and the optical lattice are successively filtered, see Fig. ||. 
For our imaging optics, the diffraction limit in resolution 
(Airy radius) is given by f.79/xm ?« 4 x A/2. Resolving 
the lattice periodicity (see later), we infer that one pixel 
of the EMCCD with a width of 16 ^m corresponds to 
294.6 nm in the object plane. 
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Figure 2: (a) Fluorescence image of atoms, (b) Vertically 
binned intensity distribution I{xi) for one region of interest, 
(c) Cumulative sum of I{xi) above the background baseline ao 
(dashed line in b). In (b), the extracted positions ^^ and the 
(scaled) fluorescence contributions ak of the atoms are shown 
(vertical lines). The shaded stripe indicates the acceptance 
region of the reliability criterion; the solid curve the resulting 
source distribution S(x) convolved with the LSF. 

The EMCCD image is taken at an exposure time of 
1 s. It provides a sampled intensity distribution lixi^ yj), 
where Xi and yj denote the horizontal and vertical po- 
sition of pixel {j,i}, respectively. The intensity distri- 
bution of a single atom trapped in the ID optical lat- 
tice shows a characteristic ellipticity, originating from the 
shape of the trapping potential, see Fig. ||(a): The atom 
is vertically confined to the lattice axis by the Gaussian 
profile of the laser beams, whereas its horizontal posi- 
tion depends on the occupied lattice site. The vertical 
width of the fluorescence spot is primarily given by the 
spread of the thcrnral wave packet of the atom in radial 
direction of the standing-wave potential. In axial direc- 
tion, the atoms are strongly confined: The horizontal 
l/y/e half width of the fluorescence spot, corresponding 



to a, 



hor 



810(±19) nm in the object plane, is mainly 
caused by diffraction within the imaging optics (with a 
theoretical value of a^°g = 647 nm). Compared with this 
width, thermal motion of the atoms (cr^h"^ = 23 nm) and 
drifts of the standing- wave (12nm/s) can be neglected. 

To simplify the axial position determination, we bin 
the intensity distribution vertically I{xi) = ^ ■ /(xi, yj), 



see Fig. ||(b). The resulting distribution is related to 
the unknown source distribution S{x) by a convolution 
equation 



HxO 



L{xi — u)S{u)du + e{xi) , 



(1) 



where L{x) is the area normalized line spread function 
(LSF) of our inraging optics and e{xi) the additive noise. 
The axial confinement of the atoms (crj^h"^ ^ "'^p'^) f^i^d 
systematical suppression of the stray light down to a ho- 
mogeneous background allow us to model S{x) as 



N 



S{x) = ao + X! "'^■'^(^ ~ ^'=) ■ 



(2) 



fc=i 



where ag denotes the constant baseline of the stray light 
background, S{x) the Dirac delta function representing 
the strongly confined atom, ak and ^k the fluorescence 
contributions and the positions of N atoms, respectively. 
Therefore, the position determination of the atoms in a 
ID lattice corresponds to a parameter estimation of the 
modeled distribution S{x). 

A sufRciently low noise is the first essential prerequi- 
site to deconvolve Eq. (P and retrieve the parameters in 
Eq. (H). In our experiment, the shot noise of the fluores- 
cence signal dominates, whereas stray light contributions 
were minimized using light traps (see Fig. nl). Readout 
noise was reduced by cooling the EMCCD. The remain- 
ing noise level is low enough to neglect its signal depen- 
dence, which greatly simplifies our nunrerical method. 

The second essential prerequisite is the precise analyt- 
ical description of the LSF. In principle, the measured 
intensity distribution of an isolated atom yields informa- 
tion on the LSF, however, a single image does not pro- 
vide the required resolution and accuracy due to noise 
and the limited EMCCD resolution. Therefore, we su- 
perinrpose the intensity distributions of up to hundred 
isolated atoms with sub-pixel accuracy and precisely fit 
the shape of the LSF. We restrict the image acquisition 
to the imaging region with negligible spatial variation of 
the LSF (« 120 ^m in object space). 

In order to extract the atomic positions in Eq. (0), we 
use a method related to the spike-convolution nrodel fit- 
ting [O which consists of three stages: segmentation, 
atom number determination and model fitting. 
In the segmentation, we divide the binned intensity dis- 
tribution into regions of interest (ROIs) which contain 
fluorescence from atoms and redundant regions which 
only contain background. From the latter we estimate 
the background baseline qq. Since parts of our numerical 
irrethod are based on Fourier transfornrs, excluding noise 
from redundant regions allows us to inrprove both the 
performance and the precision of our numerical method. 

In the following, for each ROI the same procedure 
is used. Wc determine the number of atoms N by 



cumulatively integrating the binned intensity distribu- 
tion I{xi) above the background baseline oq. Since 
each atom contributes equally, the integrated distribu- 
tion exhibits characteristic steps at integer multiples of 
the mean single atom fluorescence contribution /a, see 
Fig. 0(c). After determining the atom number from 
N = X]ieRoi(-^('^j) ~ oo)/-fa, we reestimate the (local) 
baseline ao and determine the atom positions ^k us- 
ing the trigonometric moments estimates (for details see 
Ref. [M). For each atom, the fluorescence contribution 
Ofe is determined by a least squares method, minimiz- 
ing I]ieRoi{^(^«) ~ «o - J2k=i (^kL{xi - 6)}^- Since in 
general, trigonometric moments estimates are less precise 
than maximum likelihood estimates [|l2| , we subsequently 
improve the parameters ao, a^ and ^k using them as an 
initial guess in a Levenberg-Marquardt fitting algorithm 
|_J, which then converges within few iterations. For typ- 
ical images with up to 13 atoms, we determine the atomic 
positions within less than 100 ms. 

The parameters are finally checked for passing a relia- 
bility criterion \ak — Ia\/Ia < 20 % inferred from a statis- 
tical analysis of the data jlj|. By this, we identify and 
exclude erroneous results, mainly stemming from contri- 
butions of atoms which leave the optical lattice during 
the exposure time fu^ . 

We apply our numerical method to approx. 6000 im- 
ages in order to both investigate its efficiency and accu- 
racy, and to extract the distribution of atom separations. 
For this purpose, we repeatedly load on average 4 atoms 
into the optical lattice and successively acquire several 
images of the same atom distribution. From each image 
with N atoms we determine the N < N atom positions 
in those ROIs only in which all atoms pass the reliabil- 
ity criterion and calculate the corresponding N{N — l)/2 
distances. Following Ref. Q, we additionally calculate 
the averaged distances from successively acquired images 
of the same atoms to reduce the statistical error. 

In Figure y the histograms of determined distances 
from both single images and 3 image averages are shown. 
Both reveal a periodic structure, perfectly reproducing 
the periodicity of the optical lattice, as well as a marked 
decay in the number of occurrences for small separations. 
We stress that the lattice periodicity does not enter in any 
way our estimation procedure. Thus, the strict adher- 
ence to the periodicity and the well separated histogram 
peaks provide a striking confirmation for the high preci- 
sion and reliability of our detection. For each histogram 
peak, we estimate the reliability -FngN of inferring the cor- 
rect number n of sites separating two atoms, assuming 
Gaussian distributions of the measured distances around 
d„ ~ nX/2 and fitting a sum of Gaussians centered at 
dn to the histogram. The reliability Fn is then given by 
the area of the normalized Gaussian centered at dn within 
[d„ — A/4, d„+A/4]. For site separations below the diffrac- 
tion limit (n < 5) we obtain F1-4 = 68.8-99.4%, whereas 
above the diffraction hmit ^5^29 = 97.7-99.8 %. Rcduc- 
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Figure 3: Histogram of determined distances. The upper 
rows correspond to the distances obtained from a single image; 
the lower rows to the distances obtained from the averages 
over 3 successively acquired images of the same atoms. 



ing the statistical error by averaging the distances from 3 
successively acquired images of the same atoms increases 
the rehability to F^^l > 92.0%, and F^% > 99.992%. 
This allows us to investigate possible atomic interactions 
in the nearest-neighbor regime. 

The marked decay in the number of occurrences for 
small separations can partially be traced back to the 
cases with three or more atoms occupying nearest lattice 
sites. These cases are challenging from a numerical point 
of view and cause our algorithm to provide increasingly 
inaccurate results, which then fail the reliability check. 
From simulations we deduce that this does not hold for 
pairs of nearest-neighbors separated by at least two lat- 
tice sites from other atoms. 

In order to investigate the influence of possible interac- 
tion induced effects in the physically interesting regime 
of neighboring atoms, and to get an unbiased insight in 
the statistics of atom pair separations, in the following 
experiment we focus on the distance distribution of iso- 
lated atom pairs only. For this purpose, we reduce the 
mean number of atoms loaded into the MOT to about 
two atoms using a high field-gradient of 345 G/cm, which 
also favors short distances. From single-atom images we 
infer the interaction-free position distribution, which is 
related to the overlap of the MOT and the ID lattice. 
From two-atom images, we determine the atom separa- 
tions averaging over 3 successively acquired images of 
each pair of atoms. 

In Figure 0, the histograms of single atom positions 
and pair separations are shown. The latter reveals a clear 
gap for two atoms loaded to the same lattice site, reflect- 
ing light induced atom losses. The underlying loss mech- 
anism is known | |, | |: The atoms gain sufRcient ki- 
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Figure 4: (a) Histogram of absolute positions of single atoms 
transferred from the MOT into the ID lattice. The solid curve 
shows a Gaussian fit. (b) Histogram of distances between two 
transferred atoms. Solid bars correspond to the measured dis- 
tances averaged over 3 successively acquired images; shaded 
bars to the distance distribution in terms of lattice sites. The 
solid curve shows the distribution Q{d). 



netic energy in an excited molecular potential to leave the 
trap |[rs|. For larger separations, starting from nearest- 
neighbor sites, the distance distribution follows a Gaus- 
sian shape. This shape fits well to a modeled distribution 
which assumes statistically independent positions of both 
atoms. The number of occurrences Q as a function of dis- 
tance d is then given by Q{d) = 2Qo /^ P{x)P{x+d)dx, 
where P{x) denotes the normalized fitted single atom po- 
sition distribution (Gaussian with a = 9.5 x A/2), see 
Fig. ^(a). and Qo the total number of analyzed atom 
pairs. From this fact wc conclude that, for our ID lat- 
tice, atoms in neighboring lattice sites do not affect each 
others' storage time, which is limited by background gas 
collisions. Wc stress that for isolated atom pairs with 
separations below the diffraction limit, our numerical 
procedure provides an increased reliability of more than 
95.0%, exceeding the number quoted above for clusters 
of atoms. 

Summarizing, we have presented a method to deter- 
mine the separations of individual neutral atoms in a 
ID optical lattice from a fluorescence image at all rel- 
evant distances. We have investigated the statistics of 
atom pair separations and directly observed light induced 
atom losses for atoms occupying the same lattice site. 
Our work yields promising perspectives for the imple- 
mentation and detection of controlled interactions be- 
tween atoms using, e.g., spin dependent optical potentials 



y, |2(J. Furthermore, since all stages of our numerical 
processing are extendible to higher dimensions, it points 
to high resolution imaging of 2D systems, e.g., degener- 
ate lattice gases, or nano dots on substrates in vacuum. 
For instance, finding a defect in a homogeneously filled 
Mott-insulator by inverting our method seems conceiv- 
able. 
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